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ABSTRACT 

We examine the consistency of the thermodynamics of irrotational and non-isentropic 

Tj- '. perfect fluids complying with matter conservation by looking at the integrability conditions 

of the Gibbs-Duhem relation. We show that the latter is always integrable for fluids of 

the following types: (a) static, (b) isentropic (admits a barotropic equation of state), (c) 

^ ■ the source of a spacetime for which r > 2, where r is the dimension of the orbit of the 

^ ■ isometry group. This consistency scheme is tested also in two large classes of known exact 

QQ ■ solutions for which r < 2, in general: perfect fluid Szekeres solutions (classes I and II). 

In none of these cases, the Gibbs-Duhem relation is integrable, in general, though speciflc 

K^ ' particular cases of Szekeres class II (all complying with r < 2) are identifled for which the 

pq I integrability of this relation can be achieved. We show that Szekeres class I solutions satisfy 

2 I the integrability conditions only in two trivial cases, namely the spherically symmetric 

^—1 I limiting case and the Friedman-Roberson- Walker (FRW) cosmology. Explicit forms of the 

Q\ '. state variables and equations of state linking them are given explicitly and discussed in 

O '• relation to the FRW limits of the solutions. We show that flxing free parameters in these 

I . solutions by a formal identification with FRW parameters leads, in all cases examined, 

5-1 . 
W). to unphysical temperature evolution laws, quite unrelated to those of their FRW limiting 



% 



cosmologies. 

*Work supported by CONACYT, Mexico, project No. 3567-E 



1. Introduction 

Numerous exact solutions of Einstein field equations with a perfect fluid source exist 
in the literature which do not admit a barotropic equation of state p = p{p), where p 
and p are the pressure and matter-energy density. The thermodynamic properties of these 
simple non-isentropic fluids have not been properly examined, as there is a sort of consensus 
which regards as unphysical all exact solutions with such fluid sources. However, there are 
arguments supporting the idea that barotropic equations of state are too restrictive^, and 
so it is surprising to flnd so few references in the literature^'"^'"^'^'^ studying the properties 
of these fluids. There are two related levels at which this study can be posed: (1) the 
consistency of the thermodynamic equations with the fleld equations; (2) the physical 
relevance of the state variables and equations of state linking them. Obviously, if point 
(1) is not satisfled, point (2) cannot be even addressed, also, it is quite possible to flnd 
exact solutions complying with (1) but not with (2), that is, unphysical fluids whose 
thermodynamics is formally correct. 

Coll and Ferrando"^ have formally examined and solved the question behind point (1) 
above. They derived rigurously the necessary and sufficient conditions for the integrability 
of the Gibbs-Duhem relation for perfect fluid sources complying with matter conservation. 
These conditions become a criterion to verify if a perfect fluid source of a given exact so- 
lution admits what Coll and Ferrando denote a "thermodynamic scheme" . However, these 
authors did not go beyond point (1) above (the admissibility of their consistency criterion). 
The purpose of this paper is to expand, continue and complement their work by applying 
their criterion (re-phrased in a more intuitive and practical form) to irrotational and non- 
isentropic fluids, and then to deal with point (2) mentioned above, that is, to look at the 
physical nature of the resulting equations of state in the cases when the thermodynamic 
equations are mathematically consistent. The contents of this paper are described below. 

We present in section 2 a summary of the equations of the thermodynamics of a 
general relativistic perfect fluid, together with the conditions for admissibility of a ther- 
modynamic scheme as derived by Coll and Ferrando. An immediate result is the fact that 
a thermodynamic scheme is always admissible in the following three cases: (1) the fluid is 



static; (2) the isometry group of spacetime has orbits of dimension r > 2; (3) isentropic 
fluid, admitting a barotropic equation of state p = p{p). Section 3 is concerned with the 
admissibiUty of a thermodynamic scheme for irrotational, non-isentropic perfect fluids. For 
these fluids, we investigate the conditions derived by Coll and Ferrando in terms of differ- 
ential forms expanded in a coordinate basis adapted to a comoving frame. The conditions 
of section 3 are applied in sections 4 and 5 to the irrotational perfect fluid generalization 
of Szekeres solutions (class II^~^^ and the parabolic case of class l'^'S'ii'i3,i4^^ with vanish- 
ing 4-acceleration but non-vanishing shear. These classes of solutions admit (in general) 
no isometries and are considered inhomogeneous and anisotropic generalizations of FRW 
cosmologies. An interesting result from sections 4-5 is the fact that none of these solutions 
admit a thermodynamic scheme in general, that is, with unrestricted values of the free 
parameters characterizing them. We show that for parabolic Szekeres class I solutions the 
Gibbs-Duhem relation is not integrable for r < 2. However, under suitable restrictions 
of their parameters, speciflc particular cases of Szekeres II solutions are found to admit a 
thermodynamic scheme when r < 2. For the particular cases of the solutions of sections 
4-5 found to be compatible with the thermodynamic scheme, we derive in section 6 explicit 
(though not unique) expressions of all state variables: p and p, particle number density n, 
speciflc entropy S and temperature T, as well as two-parameter equations of state linking 
them. The latter turn out to be difficult to interpret as there is no clue on how to fix 
the time dependent free parameters of the solutions. We follow a common strategy which 
consists in formally identifying these parameters with the FRW scale factor and suitable 
state variables. The resulting equations of state are totally unfamiliar and (for most cases) 
unphysical. In particular, their corresponding temperature evolution laws are unrelated to 
those of their FRW limit. Conclusions are preanted and summarized in section 7. 

We show in appendix A that the formulation of the thermodynamic scheme provided 
by Coll and Ferrando is equivalent to that given in this paper. In appendix B we show that 
parabolic Szekeres class I solutions do not admit a Killing vector of the form suggested by 
Szafron in his original paper^'^^. 

2. The thermodynamic scheme. 



Consider the perfect fluid momentum-energy tensor 



T"' = {p + p)u''u'+pg''' (1) 

where p, p and u"' are the matter-energy density, pressure and 4-velocity. This tensor 
satisfies the conservation law T"^ ;6 = which impUes the contracted Bianchi identities 

p+{p + p)e = (2a) 

h',P:b+{p + PX = (26) 

where G = «"„, iia = Ua;bU^ and h'^ = ^a ~^ UaU^ are respectively the expansion, 4- 
acceleration and projection tensor and p = u"'p,a- The thermodynamics of (1) is essentially 
contained in the matter conservation law, the condition of vanishing entropy production 
and the Gibbs-Duhem relation. The first two are given by 

(nw");a = (3a) 

{nSu'').a = (36) 

where n is the particle number density and 5" is the specific entropy. Condition (3a) 
inserted in (3b) leads to u'^S.a = 5" = 0, so that 5" is conserved along the fiuid lines. The 
Gibbs-Duhem relation can be given as the 1-form 



u) = dS — — 
T 



<^-A^ 



(4) 



where T is the temperature and d denotes the exterior derivative. The necessary and 
sufficient conditions for the integrability of (4) 

u A do; = necessary (5) 
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subjected to fulfilment of the conservation laws (2) and (3), are the conditions which Coll 
and Ferrando denote admissibility of a "thermodynamic scheme" . Using the 1-form Tiv in- 
stead of cu, these authors demonstrated that a perfect fiuid source admits a thermodynamic 
scheme iff there exists a scalar function F = F{p,p) satisfying 

e = F (6) 

the necessary and sufficient conditions for (6) to hold is, in turn, the fulfilment of the 
constraint 

{pdp - pdp) A dp A dp = (7) 

As immediate results from (6) and (7), it is easy to show that the thermodynamic 
scheme is always admissible in three important cases: (1) static fiuids; (2) isometry groups 
of spacetime have orbits of dimension r > 2 and (3) isentropic fiuids. In the first case, 
p = p = 0, and so (7) trivialy holds, while in the second case there are always local 
coordinates in which all state variables, and in particular the set {p,p,p,p) are functions 
of only two coordinates, say (t, x). Since any one of the 1-forms associated with these 
variables can be expanded as dp = p^tdt + p,xdx, etc, the wedge product in (7) vanishes. 
The third case is characterized by the existence of a barotropic equation of state p = p{p), 
implying dp = p,pdp and p = p,pP- Hence, the wedge product in (7) also vanishes. Since 
a barotropic equation of state can be shown to be equivalent to dS = and 5" = 0, the 
specific entropy in fiuids compatible with such an equation of state is a global constant 
(isentropic fiuids). For more general, non-isentropic fiuids, this quantity is a different 
constant for different observers comoving with the fiuid. Also, a relation between p and p 
(and between any other pair of state variables) necessarily involves a third state variable 
(two-parameter equations of state) . 

From these results, it is clear that perfect fiuids which might not comply with the ther- 
modynamic scheme are non-isentropic and sources of exact solutions with weaker spacetime 
symmetries: with isometry groups of dimension r < 2. Since exact solutions of this type 



are difficult to obtain, in general, we will deal in the following sections with the important 
particular case of irrotational perfect fluids. 

3. Thermodynamics of a non-isentropic irrotational perfect fluid. 

Consider an irrotational perfect fluid (vanishing vorticity tensor, iVab = 0). Since the 
4- velocity is hypersurface orthogonal, there are local comoving coordinates (t, x*), such 
that the metric, 4-velocity, 4-acceleration and projection tensor are given by 



ds^ = -N^dr + gijdx'dx^ (8a) 



u'' = N-'5t Ua = {logN),Jt i8b) 



hah = 9ijK^\ 



b 



where all the metric coefficients A^ and gij are (in general) functions of all the coordinates 
(t,a;*). In this representation, X = {l/N)X^t for all scalar functions and the Bianchi 
identities and conservation laws (2) and (3) become 

P,t + (P + P)(logv^),t = (9a) 

P,r + {p + p){logN)^, = (96) 



n = — 1=^ (9c) 



S = S{x') {9d) 

where A = det{gij) and /(a;*) appearing in (9c) is an arbitrary function denoting the 
conserved particle number distribution. Consider the coordinate basis of 1-forms (dt, dx*) 
associated with the comoving frame (8) ; the Gibbs-Duhem relation reads 
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to = S idx'^ — — 

5^ rp 



-) ^4- 

n/ ,i \n 



dx' (10) 



where the t component of u in this coordinate basis vanishes due to (9d). A sufficient 
integrabihty condition of (10) is given by 

,,, dt A dx* ,,, dx* A dx^ , , 

do; = Wu ^^ + W., ^^ = (11) 



... P[,i^,t]-n'^T[^tS,^] .fn,^Tt . nA ( Tt , n,, 

W^ti = -^ iP + P)[ — TTT--"*— - [p,z^+P,t — 

n \ n 1 n I \ 1 n 



^^ ^ P[,.n,,] + n^T[,,g,,] ^ r[,^p,j] _ (^ ^ p) /^^[,* + ^^[^ ^ ^,J] 



n T \ T J n 

where square brackets denote antisymmetrization on the corresponding indices. The nec- 
essary and sufficient condition (5) is given by 

,, dx* A dx^ A dx^ ,, dt A dx* A dx^ , , 

dc^ A u; = X,,, -^^ + Xu, -3^;^ = (12) 



Xi^ 



ijk = -P[,iP,jn^k] 

Xtij = P[,tP,in,j] 

Conditions (12) are entirely equivalent to (7) provided by Coll and Ferrando. One can 
obtain the latter form the former simply by using (2) and (3) (in their forms (9)). However, 
(11) and (12) are more intuitive than (6) and (7), as they directly incorporate state variables 
such as n, 5" and T, and their relations with p and p. Condition (12) is also more practical 
than (7), as it is easier to use n and 5" from (9c) and (9d) than to compute the set 
{p,p,p,p) in exact solutions in which these quantities can be quite cumbersome. The 
sufficient condition (11), not examined by Coll and Ferrando, is also helpful, since its 
fulfilment guarantees that (7) (or (12)) holds. As shown in the following sections, it is 



easier to test in some cases the admissibility of a thermodynamic scheme directly from 

(10) and (11). However, one must have a solution of Einstein's equations providing p and 
p in terms of the metric functions in order to verify the admissibility of the thermodynamic 
scheme. Thus we will look at this scheme in known exact solutions which are particular 
cases of (8), and to do so we suggest the following procedure: (a) solve the conditions 
(12) and substitute the solution into (10), thus identifying possible (non-unique) forms 
for S and T; insert the obtained forms of T and n into (11) in order to verify if further 
restrictions follow from the sufficient conditions. If these conditions hold, the equations 
of state linking the state variables {p,p,n, S,T) (together with their functional relation 
with respect to the metric functions) follow directly from integrating it. We shown in the 
following sections, for various known exact solutions with r < 2 and with non-isentropic 
fluid sources, that the Gibbs-Duhem relation might not be integrable and if it is integrable, 
the resulting definitions of S and T might not be unique. 

4. Geodesic perfect fluid: the Szekeres class II solutions. 

Consider the geodesic case iii = 0. In the comoving frame (8), (9) leads to A^^^ = 
and p = p{t). Without loss of generality, the metric in the comoving frame is (8a) with 
A^ = 1, so that u"" = 5 1 and all convective derivatives are simply derivatives with respect 
to the time coordinate {i.e. X = X,t for all functions X). The remaining conservation laws 
(9a), (9c) and (9d) have the same forms as given by these equations. The components of 
the 1-form associated with the Gibbs-Duhem relation (10) becomes 

while the integrability conditions da; — and uj A duj — become the particular cases of 

(11) and (12) given by 

W',. = ^(^ -pj'-) =0 (14a) 



V ^ / , V^/ , 



Wij = T[,i p,j] ^[,i^u1 = T^iS^j] = (146) 



^tij = =^ ^[,iP,j] = (15) 

The consistency between (13), its integrability conditions (14)-(15) and the field equations 
will be tested on two important classes of known exact solutions with an irrotational 
and geodesic perfect fluid source: the perfect fluid generalization of Szekeres solution^"^^. 
These solutions are divided in two subclasses: I and II. This section is devoted to the class 
II solutions, while class I solutions are investigated in the next section. 

The Szekeres class II solutions have been extensively studied as inhomogeneous gen- 
eralizations of FRW cosmologies. A usual representation of the spatial metric is given 
by 

[1 + 1(2/2 + ^2)] 

where -R(t), V{x), Vi{x), V2{x) and U{x) are arbitrary functions and k = 0, ±1. The 
function B = B{x, t) must satisfy the following constraint 

s,^p-m±p^.o (100 

which follows from the field equations. The state variables p = p{t), p and n corresponding 
to the metric (16a) are 

2RR + if + k 
p= ^2 (17a) 



_ 2RRB + 3{B + P)R'^ + k{B + 3P - V) - 2U 
^ ~ R^iB + P) ^ ^ 
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where (9c) has been used in order to arrive to (17c). The integrabihty condition (15) 
reduces to a 3-dimensional "vector product" of the type Vn x Vp — 0, whose general 
solution has the form 

p = a(i)?i^(*) + c(t) (18) 

where (a, 6, c) are arbitrary functions. It is evident that equation (18) does not hold in 
general, that is for arbitrary values of the free functions (i?, B, U, V, Vi, V2) characterizing 
the metric (16). However, it is easy to find particular cases of the latter complying with 
these equations and not admitting isometries. Given the function Bo{t) and the constant 
60, the choice of free parameters 

B = Bo{t)V, f = P, U = boV, (19a) 

leads to a solution of (15) in the form (18) with 

b{t) = 1 (196) 



R[-2RRBo + 2bo + k{l + 2Bo)] 

a{t) = (19c) 

-Dn 



where Bo{t) is restricted by (16c), and so must satisfy 

Bo + ^Bo —, -0 (20) 

However, as shown by equations (18)-(19), So must also comply with a{t) 7^ 0, or else, p 
would loose its dependence on (x, y, z) and would become the matter-energy density of a 
FRW spacetime. In order to identify T and S, we insert (19a)-(19d) into (13) yielding 



T{t) = c{t)R^Bo = [3R{R^ + k) - a{t)]Bo (21a) 
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S{x,y,z) = ^ + So (216) 

where 5*0 is an arbitrary additive constant. Regarding the sufficient condition (14): the 
part (14b) holds by virtue of T = T(t), while the part (14a) becomes (c(t) +p)T — Tp = 0, 
a condition which is identically satisfied if (20) holds, this can be verified by inserting (21a) 
and (17a) into (14a). 

Notice that satisfying the integrability conditions (14)-(15) leads to a non-unique form 
for T, since the function R is still arbitrary. Equation (20) is a second order inhomogeneous 
ordinary differential equation, having (for every choice of R) two linearly independent 
solutions. The possible forms for T, characterized by solutions of (20) , as well as equations 
of state linking p, n an p are discussed in section 6. 

5. Parabolic Szekeres class I solutions. 

A particular case of these solutions was originally examined by Szafron^^ and gener- 
alized by Bona et al.^"^, and in the framework of a two-fluid interpretation, by Sussman^^. 
The spatial metric is given in spherical coordinates x* = (r, 9, (p) as 

g.jdx'dx^ = ^^—^^i^dr^ + e^''Y^{de'^ + sin^ Odcj)^) (22o) 

where a prime denotes the derivative with respect to r and the functions Y = Y{t, r) and 
u = i^(r, 6*, (j)) are given by 

Y = {2M)^/^[v + Qw]'^/^ (226) 



€-" = !- sin^ -[i-A'^-B^- C^] + sin ^[S cos (/> + C sin (j)] (22c) 

with A(r), -B(r), C{r), M{r), and Q{r) being arbitrary functions. The pressure p = p{t) 
is determined by the otherwise arbitrary functions v — v{t) and w — w{t) as 

^v Aw 

p = --- = --- (23a) 

6V 6 w 
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while the remaining state variables p = p{t, x"^) and n = n{t, x"^) take the form 



P 



4^0 
3^0 



n 






(236) 



with 



•^ = v + wQ Q = vV + w/\ r = M ' + 3Mz/' A = TQ + 2MQ' 



(23c) 



In obtaining the expression for n in (236) we have redefined the arbitrary function /(x*) 
as/^3A//(2sin^e2'^). 

As with the class II solutions discussed previously, there are no isometries in general, 
for arbitrary values of the free parameters. This fact is proven in Appendix B, thus 
correcting the erroneous result given by Szafron^'^^, who reported in his original paper 
that a one dimensional isometry group necessarily exists in these solutions. Special cases 
of higher symmetry are obtained by specifying the free parameters, in particular a FRW 
limit follows if Y becomes separable as y = Yi{t)Y2{r). 

The necessary and sufficient condition (15) for admittance of a thermodynamic scheme,| 
with p and n given by (23b) and after lengthy algebraic manipulations, becomes (x* = 



P[,*^,. 



A{vw — wi)) 



3^303 

where 5] represents the Kronecker symbol and 



M8if,j] + ^25[>:,.] + Asz/fj,,] 



(24a) 



Ai = ^^ 



2/n',/' 



r^ + 6M^QV 



+ oog' 



(246) 



A2 = -6/M{^^[rQ' + {MQ'Y] + MQ'^(«;^)} 



(24c) 



A. 



-6M^g'^* 



(24rf) 



It follows from (24a) that the thermodynamic scheme conditions are identically satisfied if 
vw — wi) = 0. However, this relationship implies that w cc v and, according to (226), the 
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metric function Y{t, r) becomes separable. As mentioned above, this case corresponds to 
the hmiting FRW cosmology. We therefore demand that vw — wv ^ and, consequently, 
the expression in squared brackets in (24a) must vanish, i.e., 

(Ai + Aaz/")/,^ + (^2 - As/V,; = (25a) 

A3(z/:,/,^-z/,V,e) = (256) 

with T] = 6, (/). Equation (256), with A3 7^ 0, can always be satisfied by fixing the angular 
dependence of the arbitrary function /. Even if we demand the fulfillment of the simplifying 
restriction /' = 0, (256) still allows a solution which fixes the angular dependence of / and 
imposes a condition on the radial dependence of v. In fact, using (22c) and /' = it 
follows that (256) holds if B = C, B' ^ 0, and AA' = {Bq - 2B)B' where Bq is an 
arbitrary constant. Equation (25a) contains arbitrary functions of the time as well as 
spatial coordinates which can be separated by inserting (23c). Then we obtain 

vviKi + K2r^) + {vw){KiQ + K2QT + K3) + ww{KiQ^ + KsA^ + 2K3Q) = (26a) 
where 

^' = (y) ^'^'^ " 6M{/[rg' + {MQ'Y] - MQ'f'y^^ (266) 

K2 = Q'ln (26c) 

Ks = -6M^Q'^fiy[^ {26d) 

For arbitrary time functions, (26a) yields a system of three algebraic equations for the 
spatial coefficients. It then follows that this system allows a solution iff 

M'^Q'^K2 = M^Q'^f^^ = (27) 

If M = the spatial metric vanishes. Moreover, for Q' = the function Y{t, r) in (226) 
becomes separable and the spacetime is that of the FRW cosmological models. Therefore, 
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the only non trivial solution of (27) is /^^ = 0. Hence, according to (266) and (26(i), 



V , 



0, i.e. A, S, and C in (22c) must be constants and the metric (22a), after a suitable 



coordinate transformation, becomes spherically symmetric. Another possibility for solving 
(26a) is to require a linear dependence between the time functions, that is {yw) oc vi) and 
WW oc vi). However, it can easily be shown that this case leads to to oc t; and, therefore, 
we obtain the FRW limiting case. The analisis presented above shows that the Szekeres 
class I solutions admit a thermodynamic scheme in two special cases only: the spherically 
symmetric limiting case and the FRW cosmologies. 

We will now consider the spherically symmetric case which follows from (22c) through 
the conditions S = C = 0, A^ = 1, so that e" — 1. The Gibbs-Duhem relation (13) yields 



a(t) 



M' 
T 



+ W{t) 



mQ)'^ 

[ f \ 




[ f \ 



TS' 



(28a) 



where 



a{t) 






^/ N 4 9 fw 



7(t) 



4 y / w 



(286) 



y J 6 \v J 6 \w 

In order to identify the temperature and entropy from (28a), we must impose certain 
functional dependence either between the time functions or the radial coefficients. We first 
consider the former case and require that the following relationships hold: I3{t) = Kia{t) 
and 7(t) = K2Ct{t) where ki and K2 are constants. Inserting (286) into these constraints, 
we obtain a set of differential equations which may easely be integrated and yield 



w 



Kiv + eiv 



1^2 



w + e2W 



(29) 



where ei and 62 are constants of integration. Introducing the value of v into the equation 
for w and viceversa, and differentiating the resulting equations with respect to t, we get 



(k2 - 1^1) ^ - Ki{ei + e2K2)v - tit2V = 

(^2 - kI)w - Ki{ei + e2K2)w - eie2W = 
14 



(30a) 
(306) 



Since v and w are related by v/v = w/w [cf. (23a)], (30a) and (306) implies that v/v = 
w/w, a condition that leads to the FRW limiting case. Consequently, the functional 
dependence must be imposed on the radial functions contained in (28a). The constraints 



f 



Ai 



M' 

T 



and 



[ f \ 


1 

= \2 


'M'' 
[ f \ 



(31) 



with Ai and A2 being constants, can be integrated twice yielding two algebraic equations 
which imply a relationship between the arbitrary functions M and Q, namely 



Q 



X2M + T2F + (J2 



XiM + TiF + ai 
where ri, r2, o"i, and a2 are arbitrary constants and F = j f{r)dr. 

The identification of T and 5" follows by inserting (31) into (28a), thus we obtain 



(32) 



m = -3 



- +2Xiv^ - +A2«;M- 

V / \ V J \w 



M' 



(33a) 



(336) 



where Sq is an additive constant. We see from (33a) that the temperature is determined 
by the arbitrary functions v{t) and w{t) together with the arbitrary constants Ai and A2 
which must be appropriately specified for any given v{t) and w{t) in order to ensure the 
positiveness of the temperature. 

6. Equations of state and FRW limits. 

Szekeres class II solutions. 

Given a choice of R = -R(t), the constraint (20) provides -Bo(^) (or given a choice of 
the latter function, R{t) can be obtained). Since p = pit) and T = T{t), the functions 
(a(t), c(t)) appearing in (19) can always be expressed as functions of either one of the pair 
{p, T) once the constraint (20) has been solved. Thus, a generic formal equation of state 
for these solutions follows by re- writing (18) as 
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ci(t}] 
pb, n) = a{p)n + dp) = ^3(^) [i + g„(^)(g _ g„)] + <P) = P(P^ ^) (34a) 

where n has been expressed in terms of {p, S) by ehminating S from (21b) into (17c) as 

Similar forms of the equation of state, in the form p(T, n) and n{T, S) can be obtained 
by expressing p in terms of T. These formal equations of state are difficult to interpret, 
as they depend on the choice of one of the pair [R, Bq), and there is no intrinsic way to 
select these functions other than using the formal analogy between the functions R and 
p in (16a) and (17a) with the scale factor and pressure of a perfect fluid FRW cosmology 
(usually obeying a "gamma law" equation of state p = (7 — l)Ai)- This analogy implies 
fixing R{t) from solving the FRW relation 

3(7-l)(i^^ + fc) ,,fRV''' fo., 

p={'j-l)fi= — = 3(^_i)/_l (35) 

where p is given by (17a) and Rq is an integration constant. Several authors^'^"^^ have 
assumed this formal identification of these functions with their FRW analogues, hoping to 
provide a sort of "physical handle" in dealing with Szekeres class II solutions. The latter 
solutions, with these specific parameters, comply with suitably defined asymptotically 
FRW limits. However, as shown below, the thermodynamics of these fiuids bear no relation 
with that of their FRW limits, and so this is an example of how nice geometric features 
do not always correspond to nice physics. 

For Szekeres class II solutions, from equations (18) and (19), the FRW limit follows 
as a{t) -^ 0, so that p -^ pit) = c{t) = p and p{t) given by (17a) become the mass energy 
density and pressure of the limiting FRW cosmology. This FRW limit also implies S ^ Sq 
(entropy density a universal constant), and so, from (34b), we have n — > R~^ (the FRW 
form of n)). The condition a{t) -^ 0, together with (20) implies 

2RRB0 - {2kBo + k + 2bo) = 
2R'^Bo + 6RRB0 - {2kBo + k + 2bo) = 
16 



which leads to the deSitter solution if Bq 7^ 0, or FRW metrics if Bq = (for /c = ±1, 
one has Bq = —{k — 2bo)/2k, while for k = 0, Bq is an arbitrary constant, but Bq = 0). 
For the FRW limit, from (21a), the temperature function becomes T = Bq^R^, and 
H = 3{R^ + k)/R'^, so that one recovers the expected temperature law for a FRW cosmology 
with a "gamma law" equation of state, namely: T oc R^^^~^' oc n^~^''^. However, for 
Bq ^ Q and a{t) 7^ the temperature law one obtains is unrelated to these values. 
Consider first the "parabolic" case /c = 0. Inserting the FRW form 

R(t) = RQ[l + l^{t-tQ)f^'^ (36a) 

into (20) and (21a) yields the following forms for Bq and T 

D u\ I +1-97/2 260-Ro^i ('3au\ 

B„(t) = c, + c,t, - (3^ ^ 2)(37 - 2)fl§ <"'"'' 






(36c) 



where ci and C2 are integration constants of (20) and ti = 1 + |7(t — to). Notice that 
T(t) not only bears no relation with the FRW temperature law, but is a wholly unphysical 
temperature law (for 7=1, dust, T is not constant). The particular case ci = C2 = 
of (40c) was obtained by Tiomno and Lima^, who report it as unrelated to the FRW 
temperature law. However, these authors erroneously claim that the FRW temperature 
law can be recovered (for non-FRW cases) by simply redefining the equation of state 
P = p{p^'^) so that it depends on three parameters {n,p and a function of the spatial 
coordinates). Tiomno and Lima obtained T from the sufficient integrability condition 
(14a) expressed as T /T = {dp/dp)nn/n. Such integrability condition is only valid if the 
equation of state is of the two-parameter form. 

Regarding the case /c = ±1, the integration of (20) for R given by its FRW analogous 
form leads to cumbersome hypergeometric functions^^. Hence it is easier to demonstrate 
that a FRW temperature law is incompatible with such forms of R and with Bq 7^ const, 
complying with (20). Assume that T in (21a) takes the FRW form T = T^^^ = 3biR{R^ + 
k), where 61 is an arbitrary constant, equation (21a) becomes 
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2RR 2RR ^ ' 

eliminating R in (37) in terms of powers of R from (35), differentiating the result with 
respect to t and inserting the obtained forms of So and So into (20) leads, after some 
algebraic manipulation to Bq = —{k — 2bo)/2k, the relation defining the FRW limit, and 
so, indicating that a Szekeres class II solution with k = ±1 with R and T having FRW 
forms is necessarily a FRW cosmology. This means that selecting the free parameters of 
Szekeres II solutions in terms of their resemblance to FRW parameters does not lead to 
the right thermodynamics. Therefore, other criteria must be used in order to select these 
parameters. 

Szekeres class I solutions. 

As in the previous case, the physical interpretation of the thermodynamic variables 
of Szekeres' class I solutions presents serious difficulties. Consider, for instance, the tem- 
perature law (33a). Inserting the value 

of the pressure (23a) into (33a), we obtain 

T(t) = p{v^ + 2Xivw + X2W^) + -{v^ + 2Xivw + X2w'^) (38) 

This temperature law already shows an unphysical behaviour as it does not reduce, in 
general, to a constant in the limiting case of vanishing pressure. Even in the extreme case 
Ai = A2 = 0, the required behaviour holds only if i) = const. This implies an additional 
condition on the function v which cannot be satisfied in general. 

An equation of state relating p and n can be obtained from (236) by inserting the 
constraints (31) and (32), and using the definition of the entropy as in (336). Then 

4 
p = -n[{S - So){v^ + 2X1VW + X2W^) + 2tivw + Aaw^] (39) 

where the baryon number density, according to (236) and (336), can be expressed as 

n = {v + wQ)-^[{S - So){l + X2W) + r2w]-^ (40) 
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According to (38) the variable t can always be replaced by a function of one of the pair 
p and T (or both of them). Furthermore, (32) and (33) can be used to express Q as a 
function of the entropy S. Consequently, (40) may be interpreted as an equation of state of 
the form n = n{p, S) and (39) relates p with p and 5". The physical interpretation of these 
equations of state remains unclear as they explicitly depend on the choice of the arbitrary 
functions v and w. Therefore, it is necessary to fix these functions by choosing a special 
case of Szekeres class I solutions. Consider the spherically symmetric Szafron model^^'^^ 



where Rq is the constant FRW scale factor. The function M{r) remains arbitrary, and the 
thermodynamic pressure satisfies a "gamma law" and is given by 



P=(7-l)p=47nT2^ (42) 



4(7-1 ) 
3(7t)^ 

The consistency of the equations of state and the definition of temperature can be analyzed 

by looking at their specific behaviour under physically reasonable assumptions. In the FRW 

limiting case of dust (5* -^ Sq and 7 = 1), the temperature law (38) shows the correct 

behaviour as T = Iq"^ = const, while the equations of state (39) and (40) yield p = and 

n oc t~^, respectively. Obviously, this result has no relation to the evolution law expected 

at the FRW limit (n = p^^^). For radiation (7 = 4/3) the temperature evolution predicted 

by (38) becomes T oc t^/^, and the equations of state lead to p oc t~^ and n oc t~^/^ 

so that n oc p^'^; whereas the expected FRW limiting behaviour should be T oc t~^''^ 

and n oc p^'^. These examples indicate that the requirement of a thermodynamic scheme 

for the Szafron class I solutions leads to an unphysical behaviour of the thermodynamic 

variables and equations of state. 

7. Conclusions. 

We have verified the consistency of the thermodynamic equations (i.e. admittance of 
a thermodynamic scheme) for two large class of exact solutions (classes I and II Szekeres 
solutions) whose sources are non-isentropic and irrotational perfect fiuids. This work has 
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aimed at improving the study of this type of solutions, as classical fluid models generalizing 
FRW cosmologies, in contrast to a widespread attitude of simply disregarding them for 
not admiting a barotropic equation of state. 

For the particular cases of the solutions examined admitting a thermodynamic scheme, 
the resulting equations of state have an ellusive interpretation, as there is no blue print on 
how to select the free parameters of the solutions. We have show that formally identifying 
the time dependent parameter R of these solutions with the FRW scale factor leads to 
unphysical temperature evolution laws, totally unrelated to that of their limiting FRW 
cosmology. The question of how to select these parameters in a convenient matter remains 
unsolved, though the adequate theoretical framework to carry this task has been presented 
in this paper. 

It has been very interesting to find that these solutions are not compatible, in general, 
with a thermodynamic scheme. This fact seems to disqualify these solutions as classical 
fluids of physical interest. However, these exact solutions can still be useful if they are 
examined under a less restrictive framework than that of the simple perfect fluid. In this 
context, both class I and II Szekeres solutions heve been re-interpreted as mixtures of inho- 
mogeneous dust and a homogeneous perfect fluid^^'^^, subjected to adiabatic interaction. 
In all these cases, the thermodinamics is totally different (and much less restrictive) than 
that presented here. 
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Appendix A. Thermodynamic scheme conditions 

In this appendix we explicitly derive the condition for the existence of a thermody- 
namic scheme as given by Coll and Ferrando"^, and show its equivalence to the condition 

uj Adu = 0. 
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Consider a perfect fluid with mass-energy density p, pressure p and baryon number 
density n, satisfying the contracted Bianchi identities (2). If no other conditions are 
specified, (2) describe an open system. However, when an equation of state is given, the 
closure of the system is obtained. The existence of a thermodynamic scheme for the fluid 
is therefore closely related to the existence condition of an equation of state. To derive 
this condition one may use only the matter conservation law (3a). (Condition (36) cannot 
be used since it involves the entropy S which can be defined only if the fluid accepts a 
thermodynamic scheme.) 

From the matter conservation law we define the function F by 

F = e = -- (A2) 

n 



and consider the one-form 



r = dp+{p + p)dF . {A3) 



If r is completely integrable, we can define the entropy S as dS = F/T, where T is the 
integral factor which can be associated with the absolute temperature. Now we calculate 
the condition for the complete integrability of F which according to Frobenius' theorem^^ 
is equivalent to the vanishing of the three-form F A dF. Clearly, this condition will not be 
satisfied in general unless we demand certain functional dependence for F. Assuming that 
(A2) allows a solution of the form F — F{p,p), it follows that 

dF = Fpd/9 + F^pdp and F = F^p + F^pp {A4) 

where Fp = dF/dp and Fp — dF/dp. Accordingly, (A3) and the Bianchi identity (2a) 
may be written as 

F = [1 + {p + p)F^p]dp+{p + p)F,pdp (A5) 

[1 + (p + p)Fjp + (p + p)F> = (A6) 

respectively, and lead to 

p 



F = (p + p)Fp(dp-^dp) . (A7) 
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It follows from the last equation that the integrability condition for F is equivalent to 

r A dr = -(p + p)F^pp~^{pdp - pdp) A dp A dp = . (A8) 

The last expression leads to the condition (7) obtained by Coll and 

Ferrando. 

Consider now the Gibbs-Duhem relationship as given in (4), that is a; = c<;(p,p, n, T). 
Since a; is a 1-form, the condition for its complete integrability, a; A do; = 0, will be satisfied 
if w can be defined in a two-dimensional space with thermodynamic coordinates, say, p and 
p. This means that an equation of state must exist such that n = n{p,p) and T = T{p,p). 
Consequently, the Gibbs-Duhem relationship becomes 



nT 



4; Ifl - (P + P)^l dp - (p + p)'^dp] (A9) 

ni L L n J n J 



where n^p — dn/dp, etc. Introducing the Bianchi identity (2a) [with n = n{p,p)] into (A9) 

yields 

n,p fp 



uj = {p + p)^[^dp-dp\ . (AW) 

n"^! \p J 

It is now easy to see that the integrability condition uj A do; = is equivalent up to a 

constant factor to (A8). 

Appendix B. Parabolic Szekeres class I solutions and isometry groups 

In his original paper^^, Szafron introduced the particular class of parabolic Szekeres 
class I solutions that has become known in the literature as the "Szafron models". He 
claimed that the latter admit at least a one parameter isometry group, and provide specific 
forms for the components of this Killing vector. However, Szafron's claim, stated again by 
Kramer et aF, is false. We prove in this appendix that Szekeres class I solutions admit no 
Killing vector of the type suggested by Szafron. 

In section 5 we have use spherical coordinates to investigate the Szekeres class I 
solutions; however, for the investigation of their isometries it is convenient to introduce 
the original coordinates used by Szafron^ ^ and Bona et al.^"^. The spatial metric, save 
changes in notation, is given as 
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Qijdx'^dx^ = P^ 



52/3 x' 



p 



dx'^ -\ — —^{dy'^ + dz^) 



(SI) 



where a prime denotes derivative wrt a;, the function P{x, y, z) is given by 



P = W{x){y^ + z^) + Vi{x)y + V2{x) + V{x) 



(B2) 



where the functions U{x), Vi{x), V2{x) and V{x) restricted by the condition: UV — V^ 
0, while S(t, x) is determined by the field equation 



T/2 _ 1 

*^2 2 



S+fp(t)S = 



(S3) 



which can be considered as linear second order differenttial equation, and so B has the 
generic form B = M{x)v{t) +Q{x)w{t) identifying the two arbitrary integration constants 
as the arbitrary functions M{x) and Q{x). The two linearly independent solutions of (B3), 
v{t) and w{t), are related to the pressure by 



, , 4v 4w 

3v 3w 



(S4) 



Adapted to the coordinates and notation of (Bl), the generic form of the Killing vector 
provided by Szafron is 



K = y (x, y, z)dy + Z{x, y, z)dz 



(S5) 



Computing the Killing equation for the metric (Bl) and vector field (B5) leads immediately 
to the conditions 



2 r^ Zj r^ U 



(S6) 



y, + z ^ = 



,z I ^,2/ 



{Bl) 



^,z -y,y = ^ 



(BS) 
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Although Szafron's result refers to a particular case of (Bl), condition (B6) is sufficient to 
prove that this result is wrong: the components of the vector fields he provided cannot be 
components of a Killing vector of (Bl). This is so, because Szafron's particular case follows 
by giving specific forms to the functions v{t) and w{t) in (B4) and these functions are not 
involved in the calculation which leads to (B6-B9). However, even assuming Y = Y{y,z) 
and Z = Z{y,z), (B5) is not an isometry of the metric (Bl). This can be proven by 
inserting (B2) into (B9) so that the resulting equation 

(yy+.z)(^y+,(^y=o 

implies a functional dependence between the functions U and Vi of the form (Vi/P)' = 
G{y,z){U / Py , where G{y,z) is an arbitrary function of its arguments. Obviously, the 
latter condition cannot hold for an unrestricted form of P. 
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